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Abstract 

Spastic diplegic cerebral palsy(SDCP)is a common type of cerebral palsy (CP), which presents as a group of motor- 
impairment syndromes. Previous conventional MRI studies have reported abnormal structural changes in SDCP, such as 
periventricular leucomalacia. However, there are roughly 27.8% SDCP patients presenting normal appearance in 
conventional MRI, which were considered as occult SDCP. In this study, sixteen patients with occult SDCP and 16 age- 
and sex-matched healthy control subjects were collected and the data were acquired on a 3T MR system. We applied voxel- 
based morphometry (VBM) and tract-based spatial statistics (TBSS) analysis to investigate whole brain grey and white matter 
injury in occult SDCP. By using VBM method, the grey matter volume reduction was revealed in the bilateral basal ganglia 
regions, thalamus, insula, and left cerebral peduncle, whereas the white matter atrophy was found to be located in the 
posterior part of corpus callosum and right posterior corona radiata in the occult SDCP patients. By using TBSS, reduced 
fractional anisotropy (FA) values were detected in multiple white matter regions, including bilateral white matter tracts in 
prefrontal lobe, temporal lobe, internal and external capsule, corpus callosum, cingulum, thalamus, brainstem and 
cerebellum. Additionally, several regions of white matter tracts injury were found to be significantly correlated with motor 
dysfunction. These results collectively revealed the spatial patterns of whole brain grey and white matter injury in occult 
SDCP. 
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Introduction 

Cerebral palsy (CP) is a motor-impairment syndrome resulting 
from genetic or acquired disorders in the early brain development. 
About 70 to 90% of CP children presented grey and white matter 
abnormality on the conventional magnetic resonance imaging 
(MRI). Although conventional MRI is valuable in the identifica- 
tion of brain injuries, there are still 10-30% CP patients presenting 
normal appearance in conventional MRI [1,2]. We take those 
patients as occult CP. Is the brain really normal or undetectably 
abnormal in conventional MRI? 

Voxel-based morphometry (VBM) as an automated technique 
are more sensitive than conventional MRI in investigating the 
structural changes of the whole-brain [3-7] . A recent MRI study 
has employed VBM to detect regional grey matter and white 
matter volume abnormalities in CP [6]. However, previous studies 
mainly focused on CP patients who exhibited abnormal findings in 
conventional MRI. No attempt has been made to investigate grey 
and white matter changes in occult CP using this automated 



methodology. In the present study, we applied VBM to occult CP 
to explore the possible grey and white matter injury. 

With the advancement of diffusion tensor imaging (DTI) 
techniques, it has been made possible to non-invasively examine 
the whole brain white matter integrity in vivo. For the analysis of 
DTI data, fractional anisotropy (FA) is considered the most 
important diffusion metrics which reflects the degree of direction- 
ality of cellular structures within the fiber tracts by measuring 
anisotropic water diffusion [8,9]. There have been several studies 
directed at exploring the white matter integrity in CP using DTI 
method [10-12]. Unlike the previous region of interest (ROI) 
analysis method, we employed a new voxel-wise DTI analysis 
method, which is called tract-based spatial statistics (TBSS), to 
reveal the whole brain white matter injury in occult CP [13]. 

Spastic diplegic cerebral palsy(SDCP)is a common type of CP. 
Previous conventional MRI studies have reported that normal 
appearance rate of SDCP can be as high as 27.8% [2,14]. 
Therefore, we selected patients with this type of CP as the subjects. 
We hypothesized that the occult SDCP might show grey and white 
matter injury in specific regions depending on the pathology of the 
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disease. We had two objectives in mind. The first objective was to 
detect grey and white matter changes in occult SDCP using the 
whole brain VBM and TBSS technique. The second one was to 
investigate the relationship between the structural changes and 
motor dysfunction in occult SDCP. 

Materials and Methods 

Subjects 

From October 2008 to August 2013, a total of 121 SDCP 
patients aged between 2 and 14 years were included. Conventional 
MRI examinations were performed in all patients. Of the 121 
patients, 18 patients fulfilled the diagnostic criteria of occult 
SDCP. Neurological examinations were performed by two 
pediatric neurologists (DGH and ZZ) to corroborate the findings. 
The motor dysfunction scale was evaluated using the gross motor 
function classification system (GMFCS). MR images were 
interpreted and consensus about the interpretation was reached 
among 3 radiologists (MXT, WH, and MQZ). After the study was 
approved by General Hospital of Armed Police review board, 
written informed consent was obtained from all the parents 
according to the Declaration of Helsinki. 

Inclusion criteria of the occult SDCP in the study were as 
follows: 1) the age ranged from 2 to 14 years; 2) the patients must 
meet the diagnostic criteria of CP. Firstly, clinical diagnosis of a 
non-progressive motor impairment was made due to a presumably 
early insult to the developing brain; secondly, The patients may 
have one or more of the following symptoms: cognitive disabilities, 
language impairment, seizure, sensory loss, and/ or musculoskel- 
etal abnormalities; finally, objective evidence of neuromotor 
impairment, as manifested by abnormal muscle tone, strength, 
posture, reflexes, and/or motor skills, must be confirmed on 
physical examinations; 3) the patients must have spastic diplegia 
symptoms, in which the lower limbs were more severely affected 
than the upper limbs. Spasticity refers to a velocity-dependent 
increase in the muscle tone (resistance to stretch). Diplegia means 
bilateral involvement; 4) the brain showed normal-appearing in 
conventional MRI, including Tl and T2 weighted imaging [14]. 

Exclusion criteria were as follows: 1) static movement disorders 
with an onset after the age of 3, as well as progressive movement 
disorders with varying degrees of severity or psychomotor 
retardation attributable to progressive tumor, metabolic, degener- 
ative, or genetic causes; the presence of sensory deficits, ataxia, 
muscle atrophy, involuntary movements occurring in the devel- 
opment; movement disorders caused by motor neuron disease; 2) 
Data scan failure due to movement in the MRI examination. 

Eighteen patients fulfilled the above mentioned criteria and 
were included. After evaluation of the high resolution 3D 
structural images, 2 patients were excluded due to poor image 
quality. Altogether, 16 patients participated in our study and 16 
healthy children were enrolled in the normal control group. The 
inclusion criteria of the control group were as follows: full-term 
natural birth, no history of ischemia, anoxia, or genetic disorder, 
completely normal intelligence and motor function, and normal 
appearance on conventional MRI. The age and gender matched 
in the two groups (Table 1). 

Image acquisition 

A 3.0 T MR scanner (TRIO TIM, Siemens Medical Systems, 
Germany) with a SENSE 8-channel head coil was used. 
Conventional MRI study with standard imaging sequences 
preceded the VBM research protocol. 

Routine clinical pulse sequences were obtained in all subjects, 
including axial Tl -weighted fast low- angle shot sequence [repe- 



tition time (TR)/echo time (TE) = 1900/2.5 ms, slice thick- 
ness =5 mm, intersection gap = 1.5 mm, slice number = 19, field 
of view (FOV) = 220x220 mm, band width (BW) = 269 Hz, flip 
angle (FA) = 70°], sagittal and axial T2-weighted turbo spin echo 
sequence (TR/TE = 4000/87 ms, slice thickness = 5 mm, inter- 
section gap = 1.5 mm, slice number = 19, FOV = 220x220 mm, 
BW = 363 Hz, FA =140°), and axial fluid-attenuated inversion 
recovery sequence (TR/TE = 8000/ 1 1 3 ms, inversion time 
2000 ms, slice thickness = 5 mm, intersection gap = 1 .5 mm, slice 
number = 19, FOV = 220x220 mm, BW= 269 Hz, FA= 120°). 

High resolution 3D structural images were acquired by using 
magnetization prepared rapid gradient echo sequence parallel to 
the corpus callosum (TR/TE = 1900/2.58 ms, slice thick- 
ness = 0.9 mm, intersection gap = 0 mm, number of slices = 1 76, 
FOV= 240 mmx204 mm, BW= 170 Hz, FA= 15°, scan 
time = 3.38 minute, voxel size = 0.9 mmx0.9 mmx0.9 mm), cov- 
ering the entire brain. 

DTI were acquired by using spin-echo diffusion-weighted echo- 
planar imaging sequence parallel to the anterior -posterior 
commissure line (TR/TE = 4400/93 ms, slice thickness = 4 mm, 
intersection gap = 0 mm, number of slices = 32, FOV = 
220 mmx220 mm, matrix = 128x128, voxel size = 1.7 mm 
xl.7 mmx4.0 mm, MRI images were obtained from 20 non- 
collinear directions with a b value of 1000 s/mm 2 , averages times 
= 4, scan time = 6.24 minute). 

Image analysis 

The processing and data analysis were performed using the 
VBM2, an optimized VBM method [15,16], which is an extension 
toolbox of the SPM2 (Welcome Department of Cognitive 
Neurology, London, UK). A study-specific whole brain Tl- 
weighted template and prior images of grey matter, white matter, 
and cerebrospinal fluid (CSF) were created based on the Montreal 
Neurological Institute (MNI) template in SPM2. Each partici- 
pant's original image was spatially normalized based on the 
customized template and subsequendy segmented into GM, WM 
and CSF based on the customized priors. All the segmented 
images were resliced by l.Oxl.Oxl.O mm 3 voxels. Furthermore, 
voxel values in segmented images of grey and white matter were 
multiplied by the Jacobian determinants to preserve within-voxel 
volumes that may have been altered during non-linear normal- 
ization. This procedure yielded 'modulated' images, which were 
used for the group comparison of grey matter volume (GMV) and 
white matter volume (WMV). Eventually, all the GMV and WMV 
images were smoothed by a Gaussian kernel of 8 mm full width at 
half-maximum (FWHM). 

The preprocessed images were analyzed within SPM2 based on 
the framework of the general linear model [17]. In order to 
identify the difference in GMV and WMV between SDCP 
patients and the normal controls, two-sample t-test [18] was 
performed using SPM2. Brain regions with significant GMV and 
WMV changes in patients were yielded based on a voxel-level 
height threshold of p<0.001 (uncorrected) and a cluster-extent 
threshold of 50 voxels. Then the voxel-wise correlation between 
GMV, WMV and GMFCS levels was performed, with P<0.05 as 
statistically significant threshold. 

Voxel wise statistical analysis of the FA data was carried out 
using DTI-studio [19] and TBSS [13] which is a part of FSL [20]. 
FA images were created by fitting a tensor model to the raw 
diffusion data using DTI-studio. All subjects' FA data were then 
aligned into a common space using the nonlinear registration tool 
FNIRT [21,22], which utilizes a b-spline representation of the 
registration warp field [23]. The mean FA image was created and 
thinned to create a mean FA skeleton which represents the centers 
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Table 1. Demographic and clinical data of occult SDCP patients and normal controls. 





Group(n) 


Sex (M/F) 


Age(years) (mean±SD) 


Age range (years) 


Clinical factors(n) 


GMFC (n) 


Occult SDCP(16) 


9/7 


8.31 ±4.00 


2-14 


premature birth(4) 


1(8) 










low birth weight(3) 


11(4) 










perinatal asphyxia(9) 


111(4) 


NC(16) 


10/6 


7.56±3.33 


2-14 






P value 


0.4386 


0.5682 









SDCP, spastic diplegic cerebral palsy; NC, normal control; M, male; F, female; SD, standard deviation; GMFCS, gross motor function classification system. 
doi:1 0.1 371 /journal.pone.01 00451 .t001 



of all tracts common to the group. A threshold of FA > 0.2 was 
applied to exclude unmajor fiber tracts. Each subject's aligned FA 
data was then projected onto this skeleton and the resultant data 
fed into voxel wise cross-subject statistics. Statistical analysis was 
performed with TFCE method [24], which was achieved by using 
'randomise' command with '-tfce' option in FSL [25]. The 
'randomise' is a permutation test, and 5000 times permutation 
were performed to fit the group model. Regions with significant 
FA changes were yielded based on a voxel-level height threshold of 
P<0.001 (uncorrected). The ensuing voxel-wise correlation 
between FA and GMFCS levels was performed, with P<0.05 as 
statistically significant threshold. 

Results 

Demographic and clinical data of the study groups 

Descriptive characteristics and neuropsychological scores were 
shown in Table 1. No significant differences in age and gender 
were noted between the patients with occult SDCP and normal 
controls (p>0.05). Of the occult SDCP patients, 4 patients were 
born prematurely, 3 patients were with low birth weight, and 9 
patients had history of perinatal asphyxia. By GMFCS, 12 patients 
had a mild degree of motor dysfunction (Level 1 = 8, Level II = 4), 
and 4 patients had moderate degree of motor dysfunction (Level 
III = 4). 

Grey matter and white matter volume loss in occult SDCP 
using VBM method 

Grey matter volume reduction was found in the bilateral basal 
ganglia regions (lentiform nucleus and claustrum), thalamus, 
insula, and left cerebral peduncle in the occult SDCP patients, 
compared with the normal controls (Figure 1 and Table 2). In 
contrast with the normal controls, white matter volume reduction 
was located in the posterior part of corpus callosum and right 
posterior corona radiata in the occult SDCP patients (Figure 2 and 
Table 2). There was no significant grey or white matter volume 
increase in occult SDCP patients compared with the normal 
controls in this study. 

White matter injury in occult SDCP using TBSS method 

The differences of FA values between patients with occult 
SDCP and healthy controls were revealed by TBSS. Compared 
with the healthy controls, occult SDCP patients showed signifi- 
cantly reduced FA in multiple white matter regions, including 
bilateral white matter tracts in prefrontal lobe, temporal lobe, 
internal and external capsule, corpus callosum, cingulum, 
thalamus, brainstem and cerebellum (Figure 3 A-C). There were 
no significantly increased FA regions in the white matter in occult 
SDCP patients compared with the normal controls in this study. 



Correlation of the severity of white matter injury with 
GMFCS levels 

Grey and white matter volume reduction did not show 
significant correlation with motor dysfunction. FA values within 
bilateral white matter tracts of prefrontal lobe, thalamus, internal 
capsule, corpus callosum and brainstem showed a negative 
correlation with GMFCS levels (P<0.05) (Figure 4) 

Discussion 

Patients with SDCP whose brain abnormalities are unobserv- 
able on conventional MRI are termed "occult SDCP" in this study 
[1]. However, subde structural injury may exist in these patients, 
which, supposedly, is closely related to motor dysfunction. We 
combined two highly sensitive techniques including VBM and 
TBSS to reveal the grey and white matter abnormalities in occult 
SDCP in the hope of offering a clearer understanding of the 
disease. By using VBM method, the occult SDCP showed grey 
matter reduction in multiple brain regions mostly concentrated in 
basal ganglia, thalamus, insula, while white matter reduction 
tended to be found in posterior corpus callosum and the right 
corona radiata. As shown by TBSS, occult SDCP exhibited 
significandy reduced FA values in several white matter tracts, some 
of which including prefrontal lobe, thalamus, internal capsule, 
corpus callosum and brainstem showed significant negative 
correlation with motor dysfunction scores. 

Previous studies have shown grey matter volume reductions in 
the CP patients, which suggested neuronal degeneration and 




Figure 1 . Grey matter volume differences between occult SDCP 
patients and normal controls by VBM analysis. Occult SDCP 
patients showed significantly decreased grey matter volume in bilateral 
lentiform nucleus, claustrum, thalamus, insula, and left cerebral 
peduncle. Reader's left is subject's left. 
doi:1 0.1 371/journal.pone.01 00451 .g001 
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Table 2. Areas of grey and white matter volume loss in occult SDCP patients compared with NC. 





Anatomic area 


Talairach coordinate 




Cluster size 


Z scores 


P value 




X 


y 


z 








Grey matter 


L cerebral peduncle 


-7 


-15 


-21 


136 


3.67 


0.000 


R-L thalamus 


-13 


-25 


19 


6437 


3.95 


0.000 




18 


-25 


20 




3.91 


0.000 




-2 


-26 


11 




3.73 


0.000 


R lentiform nucleus/claustrum/insula 


31 


-14 


10 


1538 


3.48 


0.000 




30 


2 


13 




3.37 


0.000 


L lentiform nucleus 


-29 


-14 


9 


302 


3.34 


0.000 


L lentiform nucleus/claustrum/insula 


-26 


8 


10 


765 


3.39 


0.000 


White matter 


R posterior corona radiata 


25 


-32 


21 


1167 


3.47 


0.000 


corpus callosum 


4 


-23 


23 


410 


3.4 


0.000 




0 


-28 


18 




3.2 


0.001 



NC, normal control; R, right; L, left. 
doi:1 0.1 371 /journal.pone.01 00451 .t002 



Figure 2. White matter volume differences between occult 
SDCP patients and normal controls by VBM analysis. Occult 
SDCP patients showed significantly decreased white matter volume in 
posterior part of corpus callosum and right posterior corona radiata. 
Reader's left is subject's left. 
doi:1 0.1 371 /journal.pone.01 00451 .g002 



damage [26-31], As the basal ganglia regions are very sensitive to 
the prenatal hypoxic-ischemic insult, these regions are vulnerable 
to injuries and degeneration, as was confirmed by recent studies 
[6,1 1]. The thalamus atrophy is also an important finding in occult 
SDCP patients. Recent neuropathological studies indicated the 
damage of thalamus in CP due to diffuse neuronal loss and gliosis 
[32]. In previous CP studies, several quantitative volumetric MRI 
studies demonstrated the decreased volume of thalamus in CP 
patients [27,33,34]. Furthermore, DTI measurements of the 
thalamus also revealed a significant reduction of FA value in 
bilateral thalamus [10,1 1]. In addition, insula was also found to be 
involved. As insula is adjacent to the basal ganglia regions where 
connections exist, its volume reduction might result from neuronal 
degeneration occasioned by hypoxic-ischemic insult. Nevertheless, 
we did not find grey matter atrophy in other regions, such as 
sensorimotor cortex, frontal, parietal, or occipital lobes, which 
finding was different from those reported by previous studies. 
Considering that only occult SDCP was studied, we speculated 
that the basal ganglia, thalamus and insula were preferential and 
sensitive regions, which, presumably, are the very early affected 
sites in SDCP. 

By using VBM method, we also found white matter injuries in 
occult SDCP, which were chiefly located in the posterior corpus 




Figure 3. TBSS results of FA values differences between occult SDCP patients and normal controls groups. Green is the mean white 
matter skeleton of all subjects; red is the regions with reduced FA in occult SDCP patients. TBSS results for increased FA in occult SDCP patients did 
not show any statistically significant differences at P<0.001 (uncorrected). (A) axial position; (B) sagittal position; (C) coronal position. Reader's left is 
subject's left. 

doi:1 0.1 371 /journal.pone.01 00451 .g003 
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Figure 4. Voxel-wise correlation analysis between FA values 
and GMFCS levels on TBSS. Green is the mean white matter 
skeleton of all subjects; Areas in blue represent regions showing 
significant negative correlation between FA values and GMFCS levels. 
Reader's left is subject's left. 
doi:1 0.1 371 /journal.pone.01 00451 .g004 

callosum and right posterior corona radiata. Previous studies 
demonstrated that white matter pathology was an early marker of 
the disease, which mainly presented with gliosis and microscopic 
cystic changes resulting from hypoxic-ischemic insult [11]. 
According to previous CP studies, diminished motor cortical 
connectivity and the corticospinal tract could be the pathophys- 
iological mechanism responsible for motor dysfunction. Because 
the bilateral motor cortex is connected by the posterior part of 
corpus callosam, and posterior corona radiata is part of parieto- 
occipital white matter tract, we speculated that posterior corpus 
callosum and corona radiata might be the predilection white 
matter injury site of SDCP. 

One of the advantages of DTI over VBM, which is capable of 
identifying the subtle white matter injury, is that it can provide 
fuller and more detailed information on white matter microstruc- 
ture injury based on the diffusion properties. In this study, we 
compared the group differences in FA values between occult 
SDCP patients and healthy controls using a novel TBSS method, 
and found the diffuse white matter injury in multiple brain regions 
such as prefrontal lobe, temporal lobe, internal and external 
capsule, corpus callosum, cingulum, thalamus, brainstem and 
cerebellum. A previous DTI study showed marked reduction of 
white matter in CP patients mainly in the internal capsule, 
posterior thalamic radiation, superior corona radiata and com- 
missural fibers [35] . Another study using DTI method showed FA 
values were significantly lower in posterior white matter tracts and 
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